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ABSTRACT 

We present a new model for computing consistently the line and continuum emission 
from galaxies, based on a combination of recent population synthesis and photoion- 
ization codes. We use effective parameters to describe the H n regions and the diffuse 
gas ionized by single stellar generations in a galaxy, among which the most important 
ones are the zero-age effective ionization parameter, the effective gas metallicity, and 
the effective dust-to-heavy element ratio. We calibrate the nebular properties of our 
model using the observed [Oiii]/H/3, [Oii]/[Oiii], [Sn]/Ha, and [Nll]/[Sn] ratios of 
a representative sample of nearby spiral and irregular, starburst, and Hn galaxies. 
To compute whole (line plus continuum) spectral energy distributions, we include the 
absorption by dust in the neutral interstellar medium (ISM) using a recent simple 
prescription, which is consistent with observations of nearby starburst galaxies. Our 
model enables us to interpret quantitatively the observed optical spectra of galaxies in 
terms of stars, gas, and dust parameters. We find that the range of ionized-gas prop- 
erties spanned by nearby galaxies implies factors of 3.5 and 14 variations in the Ha 
and [On] luminosities produced per unit star formation rate (SFR). When account- 
ing for stellar Ha absorption and absorption by dust in the neutral ISM, the actual 
uncertainties in SFR estimates based on the emergent Ha and [O n] luminosities are 
as high as several decades. We derive new estimators of the SFR, the gas-phase oxy- 
gen abundance, and the effective absorption optical depth of the dust in galaxies. We 
show that, with the help of other lines such as [On], H/3, [Om] [Nil], or [Sn], the 
uncertainties in SFR estimates based on Ha can be reduced to a factor of only 2-3, 
even if the Ha line is blended with the adjacent [Nil] lines. Without Ha, however, 
the SFR is difficult to estimate from the [Oil], H/3, and [Om] lines. The reason for 
this is that the absorption by dust in the neutral ISM and the ionized-gas parameters 
are then difficult to constrain independently. This suggests that, while insufficient by 
itself, the Ha line is essential for estimating the star formation rate from the optical 
emission of a galaxy. 

Key words: galaxies: general - galaxies: ISM - galaxies: stellar content. 



1 INTRODUCTION 

The most prominent tracers of star formation in the inte- 
grated spectra of galaxies are the emission lines produced 
by short-lived, massive stars. Several approximate formulae 
have been derived to convert the Ha and [O n] luminosities 
of observed galaxies into star formation rates (SFRs; see the 
review by Kennicutt 1998). Calibrations of the Ha luminos- 
ity emitted per unit SFR have consisted usually in applying 
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dust-free case B recombination to the stellar ionizing radia- 
tion predicted by a population synthesis model (e.g., Kenni- 
cutt 1983). The scaling of this luminosity by some empirical 
[On] /Ha ratio provides a calibration of the [On] luminos- 
ity emitted per unit SFR (Gallagher, Hunter & Bushouse 
1989; Kennicutt 1992b; see Barbaro & Poggianti 1997 for a 
more refined calibration) . In reality, the absorption of ioniz- 
ing photons by dust in H II regions could affect seriously the 
Ha luminosity emerging from a star- forming galaxy (e.g., 
Mathis 1986b). Moreover, since the [On]/Ha ratio can vary 
by an order of magnitude in nearby star-forming galaxies, 
SFR estimates based on this ratio are subject to substan- 
tial uncertainties (Kennicutt 1998). To address these issues, 
we require a model for relating consistently the luminosities 
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of emission lines to the stars, gas, and dust parameters of 
galaxies. 

The optical-line ratios in the integrated spectra of 
nearby star-forming galaxies are similar to those in the spec- 
tra of individual H II regions (e.g., Lehnert & Heckman 1994; 
Kobulnicky, Kennicutt & Pizagno 1999). Photoionization 
models have proved essential to interpret the emission-line 
properties of H II regions in terms of stars and gas param- 
eters (e.g., Evans & Dopita 1985). Thus, these models can 
potentially help us also interpret the nebular emission from 
galaxies in terms of macroscopic star formation parameters. 
In fact, photoionization models have already been combined 
with evolutionary population synthesis codes to constrain 
the ages and metallicities of giant extragalactic Hn regions 
and Hn galaxies (e.g., Garcia- Vargas, Bressan & Diaz 1995, 
Stasiriska & Leitherer 1996, and references therein). These 
studies, however, are not suited to spectral analyses of galax- 
ies containing several stellar generations. Furthermore, un- 
til recently, combinations of photoionization and popula- 
tion synthesis models were limited by the lack of a simple 
prescription to compute consistently the effects of dust on 
the line and continuum emission from galaxies (see however 
Chariot & Fall 2000). It is now possible to build a model 
to interpret simultaneously the signatures of stars, gas, and 
dust in the integrated spectra of galaxies. The purpose of 
this paper is to present such a model. 

We combine recent population synthesis and photoion- 
ization codes to compute the line and continuum emission 
from galaxies. We use effective parameters to describe the 
H II regions and the diffuse gas ionized by single stellar gen- 
erations in a galaxy. We first calibrate the nebular proper- 
ties of our model using the observed [O iii]/H/3, [O n]/[0 ill], 
[S ll]/Ha, and [N ll]/[S n] ratios of a representative sample of 
nearby spiral and irregular, starburst, and Hn galaxies. To 
compute whole (line plus continuum) spectral energy dis- 
tributions, we then include the absorption by dust in the 
neutral interstellar medium (ISM) using a simple prescrip- 
tion, which is consistent with observations of nearby star- 
burst galaxies (Chariot & Fall 2000) . Our model succeeds in 
accounting quantitatively for the optical line and continuum 
emission from nearby star- forming galaxies of various types. 
Such spectral fits provide insights simultaneously into the 
star formation history, metallicity, and absorption by dust 
in the galaxies. This enables us to identify and calibrate new 
estimators of the star formation rate, the gas-phase oxygen 
abundance, and the effective absorption optical depth of the 
dust in galaxies. 

We present our model in §2, where we express the nebu- 
lar emission of a galaxy in terms of effective gas parameters. 
In §3, we compare our model with observations and identify 
the specific influence of each parameter on the integrated 
spectral properties of galaxies. In §4, we use the observed 
relations between various integrated spectral properties of 
nearby star-forming galaxies to construct estimators of the 
star formation rate, the gas-phase oxygen abundance, and 
the effective absorption optical depth of the dust. We express 
our results in terms of simple formulae derived for different 
assumptions about the available spectral information. Our 
conclusions are summarized in §5. 



2 THE MODEL 

In this section, we describe our model for the emission from 
stars and gas in galaxies. To begin with, in §2.1-2.3, we fo- 
cus on the production of stellar radiation and its conversion 
into emission lines and recombination continuum by the gas. 
We compute the luminosity emerging from a galaxy in §2.4, 
where we include the absorption of photons emanating from 
the ionized regions by dust in the neutral (and molecular) 
ISM. We do not consider here the contribution to the ioniz- 
ing radiation by active galactic nuclei (AGNs; see §4). 

2.1 Basic Assumptions 

The luminosity per unit wavelength (t) produced at time 
t by the stars and the ionized gas in a galaxy can be ex- 
pressed generally as 

L+(t)= [ dt'i>(t-t')Sx(t')T+(t,t'). (1) 
Jo 

Here ip(t — t') is the star formation rate at time t — t' , S\(t') 
is the luminosity emitted per unit wavelength and per unit 
mass by a stellar generation of age t' , and T^{t,t') is the 
'transmission function' of the ionized gas, defined as the frac- 
tion of the radiation produced at wavelength A at time t by 
a generation of stars of age t' that is transferred by the 
gas these stars ionize. Thus, T^(t,t') must be regarded as 
the average transmission of the gas ionized by all stars of 
age t' within the galaxy. If the ionized regions are bounded 
by neutral material, T+(t,t') will be close to zero at wave- 
lengths blueward of the H-Lyman limit but greater than 
unity at the wavelengths corresponding to emission lines. 
To compute L~^(t) from equation ([[]), we make the following 
simplifying assumptions. 

We neglect the ionization of the gas by hot 
intermediate- and low-mass stars; we thus write 

T+(t,t') = l, fort'>i OB , (2) 

where t OB is the lifetime of massive stars producing most 
of the ionizing photons (§2.2). At ages greater than t OB , 
planetary nebula nuclei are the dominant sources of ionizing 
photons in a single stellar generation. The ionizing radiation 
from planetary nebula nuclei, however, is typically less than 
0.1 per cent of that produced at earlier ages by massive stars, 
even for initial mass functions (IMFs) with steep slopes and 
low upper cutoff masses (e.g., Chariot & Fall 1993; Binette 
et al. 1994). Since we are interested in galaxies in the process 
of forming massive stars, we therefore neglect the influence 
of the radiation from old planetary nebula nuclei on the gas. 

We also neglect the influence of stellar winds and super- 
novae on the emission from galaxies. This is justified by both 
models and observations of star-forming galaxies. Leitherer 
& Heckman (1995) have computed the power injected into 
the ISM by stellar winds and supernovae for stellar popu- 
lations with various IMFs and star formation histories. For 
a normal IMF, this amounts to less than 10 per cent of 
the power of the ionizing radiation from massive stars (see 
Fig. 69 of Leitherer & Heckman 1995). Observations sup- 
port these results. About 20 to 50 per cent of the H-Balmer 
line emission from nearby spiral and irregular galaxies arises 
from 'diffuse ionized gas', which is not directly associated 
to individual Hn regions (e.g., Hunter & Gallagher 1990; 
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Hoopes, Walterbos & Greenwalt 1996; Ferguson et al. 1996; 
Martin 1997). This gas, however, appears to also be heated 
primarily by the ionizing radiation from massive stars rather 
than by shocks (Mathis 1986a; Hunter & Gallagher 1990; 
Hoopes et al. 1996; Martin 1997; Wang, Heckman & Lehnert 
1997). Direct evidence that 20 to 50 per cent of the ioniz- 
ing photons produced in individual Hn regions might leak 
into the ambient ISM further supports this result (Oey & 
Kennicutt 1997). Hence, for most purposes, the influence of 
stellar winds and supernovae on the nebular emission from 
galaxies can probably be neglected. 

Finally, we assume for simplicity that the transmission 
function of the ionized gas depends only on the age t' of the 
stars producing the ionizing radiation, i.e. 

T+(t,t')=T+(t'), fort'<t OB . (3) 

This approximation is adequate for the purpose of interpret- 
ing observations of star-forming galaxies. In fact, since the 
nebular emission depends only on the current gas content of 
the galaxies, specifying the dependence of on the galaxy 
age t would involve arbitrary assumptions about the chemi- 
cal enrichment history. We now appeal to standard codes to 
compute the radiation S\(t') from the stars and the trans- 
mission function T^(t') of the gas in our model. 



2.2 Stellar Radiation 

We compute the luminosity S\(t') emitted per unit wave- 
length and per unit mass by a stellar generation of age t' 
using the most recent version of the Bruzual & Chariot 
(1993) population synthesis code. This includes all phases 
of stellar evolution, from the zero-age main sequence to su- 
pernova explosions for progenitors more massive than 8 Mq , 
and to the end of the white dwarf cooling sequence for less 
massive progenitors. Models can be computed at metallic- 
ities Z = 0.0001, 0.0004, 0.004, 0.008, 0.02, 0.05, and 0.1 
(Zq = 0.02). In the models used here, stars evolve along 
the Padova tracks (Alongi et al. 1993; Bressan et al. 1993; 
Fagotto et al. 1994a,b,c; Girardi et al. 1996). The tracks in- 
clude mild overshooting in the convective cores of stars more 
massive than 1.0 Mq (with a reduced efficiency between 1.0 
and 1.5 Mq), as indicated by various observations of Galac- 
tic star clusters. We adopt the semi-empirical library of stel- 
lar spectra compiled by Lejeune et al. (1997, 1998) to de- 
scribe the emission from stars of all metallicities. This library 
relies on model atmospheres by Kurucz (1992, 1995) for the 
hotter (O-K) stars with effective temperatures in the range 
3, 500 < T c ff < 50, 000 K, complemented with models by 
Bessell et al. (1989, 1991) and Fluks et al. (1994) for cooler 
M giants, and with models by Allard & Hauschildt (1995) for 
M dwarfs. Lejeune et al. (1997, 1998) used empirical correc- 
tions to refine all these models at optical and infrared wave- 
lengths. The spectra of stars hotter than T e g = 50, 000 K 
in the population synthesis models (i.e., a few short-lived 
Wolf-Rayet stars and central stars of planetary nebulae) 
are approximated by pseudo-blackbodies, which include the 
discontinuities across the important photoionization edges 
blueward of the H- Lyman limit (as parameterized by Shields 
& Searle 1978 to reproduce the models of Hummer & Mi- 
halas 1970). Nearly all the ionizing radiation from a stellar 
population with a normal IMF is produced by stars with 



Tea <; 30, 000 K. Thus, the results presented in this paper 
rely heavily on the Kurucz model atmospheres. 

The Kurucz (1992, 1995) models are static, plane- 
parallel atmospheres in local thermodynamic equilibrium 
(LTE) with line blanketing by over 58 million atomic and 
molecular transitions. They have been tested against other 
recent models that include non-LTE effects, spherical geome- 
try, atmospheric expansion, and stellar winds, but which are 
generally restricted to smaller ranges of temperature and in- 
clude line blanketing by fewer transitions (Vacca, Garmany 
& Shull 1996, and references therein; Schaerer & de Koter 
1997, and references therein). These comparisons indicate 
that the Kurucz models produce more H-ionizing photons 
(A < 912 A) than spherical models including non-LTE and 
wind effects, by an amount that increases from a few percent 
at T eS = 50, 000 K to about 25 per cent at T eS = 35, 000 K 
for main-sequence stars. The Kurucz models produce also 
less Hel-ionizing photons (A < 504 A), by an amount that 
decreases from 40 to 25 per cent in this temperature range. 
Thus, for a stellar population with a normal IMF, we expect 
the overall discrepancy in the H- and Hel-ionizing continua 
produced by the two types of models to be less than about 
15 and 30 per cent, respectively. This level of accuracy is 
sufficient for our purposes. In fact, the absolute accuracy 
of the Kurucz models should be higher, since the effects 
of line blanketing are probably underestimated in spheri- 
cal models including non-LTE and wind effects. Unfortu- 
nately, observations of the H-ionizing continuum of stars 
with T e ff ^ 30, 000 K are not available, while the spectra 
of two stars with T e fr <; 20, 000 K have led to ambiguous 
interpretations (Cassinelli et al. 1995; Vacca et al. 1996; 
Schaerer & de Koter 1997). We note that the neglect of wind 
effects in the Kurucz models may result in a deficiency of 
He n-ionizing photons (A < 228 A; see Schaerer & de Koter 
1997). Therefore, we do not rely on the predictions of the 
population synthesis code at wavelengths blueward of the 
He II- ionizing edge. 

The main adjustable parameters in the population syn- 
thesis code are the IMF, the star formation rate, and the 
metallicity. For simplicity, we assume that the IMF can be 
approximated by a single power law, 4>( m ) rn~ 1 ~ x [defined 
such that (f>(m)dm is the number of stars born with masses 
between m and m + dm] , with an exponent x and lower and 
upper cutoff masses mi and ran. We include mo as an ad- 
justable parameter and, unless otherwise indicated, adopt a 
Salpeter IMF (a; = 1.35) with ra L = 0.1 Mq . We compute 
a model spectrum at ages from 1 x 10 5 yr to 2 x 10 10 yr in 
time steps that increase from 1.5 x 10 4 yr to 2.5 x 10 8 yr at 
late ages. At each time step, the spectrum is defined over the 
range of wavelengths from 91 A to 160 /im, with samplings 
of 10 A in the ultraviolet, 20 A in the visual, and increasing 
further in the infrared. The lifetimes of the stars producing 
most of the ionizing radiation are typically 3-5 xlO 6 yr at all 
metallicities. To be conservative, we adopt i OB = 1 x 10 7 yr 
in equations (^) and (^), corresponding to a drop by over 99 
per cent in the production rate of ionizing photons. 

We are also interested in the absorption equivalent 
widths of the low-order Balmer lines of Hi. The spectral 
resolution of the population synthesis code is too low for 
line strengths to be measured reliably. We adopt the stan- 
dard procedure that consists in parameterizing absorption- 
line strengths as functions of the stellar effective temper- 
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ature, gravity, and metallicity (e.g., Worthey et al. 1994). 
This method has been developed mainly to interpret the 
features of late-type dwarf and giant stars in the spectra of 
old stellar populations. As a consequence, appropriate cali- 
brations do not exist for the Ha and H/3 absorption equiva- 
lent widths of hot-main sequence and supergiant stars that 
dominate the spectra of young stellar populations. Cananzi, 
Augarde & Lequeux (1993) did collect accurate H7 and H<5 
equivalent widths of nearby stars over complete ranges of 
spectral types and luminosity classes for the specific pur- 
pose of population synthesis analyses. Detailed Balmer-line 
models by Kurucz (1992) indicate that the Ha, H/3, and H7 
equivalent widths of a star are always the same to within 
30 per cent and do not depend sensitively on metallicity 
(see also Gonzalez-Delgado & Leitherer 1999) . We therefore 
adopt the H7 absorption equivalent widths listed in Table 2 
of Cananzi et al. (1993) as estimates of the Ha and H/3 
equivalent widths of stars of all metallicities. We compute 
the line equivalent widths in the integrated spectra of stellar 
populations by weighting the contributions from individual 
stars by their level of continuum (see, e.g., Bressan, Chiosi & 
Tantalo 1996 for a description of this standard procedure). 

2.3 Transmission by the Photoionized Gas 

We compute the transmission function T^(t') of the gas ion- 
ized by stars of age t' (eqs. Q]-||) using the standard pho- 
toionization code CLOUDY (version C90.04; Ferland 1996). 
This requires us to specify the parameters of the gas. Here, 
we adopt 'effective' parameters to describe the ensemble of 
Hn regions and the diffuse gas ionized by all stars of age 
t' within a galaxy. This is motivated by the fact that the 
optical-line ratios in the integrated spectra of nearby spi- 
ral galaxies are similar to those in the spectra of individ- 
ual Hn regions (Kobulnicky et al. 1999; see also §3 below). 
Thus, photoionization models that reproduce the properties 
of H 11 regions should be adequate, modulo some adjustment 
of the parameters, to reproduce the integrated properties of 
whole galaxies. For example, at fixed metallicity, the ratio 
of ionizing-photon to gas densities is expected to be lower 
when averaged over a whole galaxy than in individual H 11 
regions because of the contribution by diffuse ionized gas 
(e.g., Lehnert & Heckman 1994; Martin 1997). 

In the photoionization code CLOUDY, the gas is de- 
scribed as spherical concentric layers centred on the ioniz- 
ing source (assumed to be pointlike). The ratio between the 
radius of the innermost layer r ln and the Stromgren radius 
Rs fixes the actual geometry of a model and the ionization 
profile of the gas. The Stromgren radius is defined by 



R 3 S = 3Q/(4nn^ea B ) 



(4) 



where Q is the rate of ionizing photons, 71h the hydrogen 
density, e the volume-filling factor of the gas (i.e., the ratio 
of the volume-averaged hydrogen density to nu), and qb 
the case-B hydrogen recombination coefficient (Osterbrock 
1989). The 'ionization parameter' is defined as the ratio of 
ionizing-photon to gas densities at the distance r from the 
ionizing source, 



U(r) = Q/(47rr n H c) , 



(5) 



where c is the speed of light. For n n <; _Rs , the total thickness 
Ar of the ionized region is much smaller than the Stromgren 



radius, Ar -C Rs, and the ionization parameter is roughly 
constant throughout this region, U(r) « U(ri n ). The geom- 
etry in this case is approximately plane-parallel, and the 
volume-averaged ionization parameter is (U) ~ U(ri n ). For 
Tin <C Rs, however, U(r) is a strong function of r, and the 
thickness of the ionized region approaches the Stromgren ra- 
dius, Ar ~ Rs- The geometry in this case is truly spherical, 
and the volume- averaged ionization parameter is, assuming 
constant %, 



(U) « 3Q/(47r#in H c) = W(R S ) 



(6) 



Substituting Rs from equation (^) and neglecting the weak 
dependence of ob on r through the electron temperature, 
this yields 



\ 47T 



We find that, at fixed (U), our results are not sensitive to the 
assumed geometry of the gas. For simplicity, we adopt ri n = 
0.01 pc and hence spherical geometry in all the following 
applications. 

In our approach, the parameters of the photoionization 
code are effective ones, which describe the ensemble of H 11 
regions and the diffuse gas ionized by a single stellar genera- 
tion in a galaxy. Since the ionizing radiation changes as the 
stars evolve, these effective parameters depend on time. We 
call Q(t') the effective rate of ionizing photons seen by the 
gas irradiated by stars of age t' throughout the galaxy. For 
a fixed IMF, Q(t') defines a characteristic mass M* through 
the relation 



Q(t') 



he 



d\\S x (t') 



(8) 



where S\(t') is again the luminosity emitted per unit wave- 
length and per unit mass by a stellar generation of age t 
(\l = 912 A). The quantity M* should be interpreted as 
the effective mass of the ionizing star clusters in the galaxy. 
For simplicity, we assume that the effective density % and 
the effective volume-filling factor e of the ionized gas do not 
depend on t' . By analogy with equations (^)-(0), we define 
the effective ionization parameter of the gas irradiated by 
stars of age t' and its average over volume as 



U(t',r) = Q{t')/{4mr 2 nnc), 
2/3 

(U)(t')^3U(t',R s ) « 

c 



3 Q(t r 



e 2 n H 



4tt 



1/3 



(9) 



(10) 



Here r and Tis pertain to the Hn region ionized by a star 
cluster of effective mass M*. In our model, therefore, a 
galaxy containing several stellar generations comprises dif- 
ferent gas components having different effective ionization 
parameters. 

We must also specify the abundances of heavy ele- 
ments and dust in the gas. The 'cosmic' abundances of 
heavy elements at solar metallicity are constrained to within 
only a factor of two by observations of Galactic stars (see 
the discussions by Garnett et al. 1995 and Snow & Witt 
1996). As a result, various assumptions can be found in 
previous photoionization studies about the cosmic abun- 
dances of important gas 'coolants' such as N, O, and S 
(e.g., McCall, Rybski & Shields 1985; Campbell 1988; Mc- 
Gaugh 1991; Shields & Kennicutt 1995; Stasinska & Lei- 
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therer 1996; Tresse et al. 1996; Martin 1997; Bresolin, Ken- 
nicutt & Garnett 1999). Here, we adopt the following cosmic 
abundances for these elements by number relative to hydro- 
gen: (N/H) cosmic = 5.0xl0" 5 , (0/H) cosmic = 6.6xl0" 4 , and 
(S/H) cosm i c = 1.3 x 10~ 5 . These choices, which are within 
the ranges considered in the above studies, were motivated 
in part by the desire to match roughly the observed typical 
(i.e., median) properties of the samples of nearby galaxies 
and Hn regions investigated in §3. For all the other elements, 
we adopt the solar-composition abundances of Grevesse & 
Anders (1989; as extended by Grevesse & Noels 1993). Fur- 
thermore, for simplicity, we assume that the abundances of 
all elements except nitrogen scale linearly with the effec- 
tive gas metallicity, noted Z. Observations of Hn regions in 
the Milky Way and other nearby galaxies indicate that, at 
high metallicity, nitrogen is enhanced by secondary nucle- 
osynthetic processing (e.g., Vila-Costas & Edmunds 1993; 
Henry & Worthey 1999, and references therein). We assume 
that N/H scales as Z 2 for Z > OAZq, consistent with the 
observations compiled in Fig. 6 of Henry & Worthey (1999). 
In all applications in this paper, we adopt the same metal- 
licity for the gas as for the stars (in practice, we adopt the 
available stellar metallicity closest to that set for the gas). 

In our model, we treat the mass fraction of heavy el- 
ements locked into dust grains in the ionized gas as a free 
parameter. We refer to this as the 'effective dust-to-heavy 
element ratio', noted £d- For simplicity, we assume that, at 
fixed metallicity, the fraction of each refractory element de- 
pleted onto dust grains from the gas phase scales linearly 
with £d (the secondary production of non-refractory nitrogen 
implies that, when Z changes, this fraction scales almost but 
not perfectly linearly with £d). For all the elements except 
silicon, we adopt the default depletion factors of CLOUDY 
for the average ISM (Cowie & Songaila 1986; Jenkins 1987). 
These factors correspond to a default dust-to-heavy element 
ratio of 0.46 for the above cosmic abundances. We adopt 
the more recent estimate of the silicon depletion factor from 
Hubble Space Telescope (HST) observations of H II regions by 
Garnett et al. (1995). We take this to correspond to a dust- 
to-heavy element ratio of 0.30, which is expected to be more 
typical of Hn regions, where grains can be eroded and de- 
stroyed (see Garnett et al. 1995; Shields & Kennicutt 1995). 
The optical properties of the dust in the photoionization 
code are based on the standard Draine & Lee (1984) grain 
model with extensions at ionizing wavelengths by Martin 
& Rouleau (1989). The resulting extinction curve displays 
a maximum near 17 eV (A pb 730 A). The most important 
effect of dust in the photoionized gas is the depletion of 
coolants (Shields & Kennicutt 1995). Other effects include 
the absorption and scattering by dust grains of the incident 
radiation and the photoelectric heating and collisional cool- 
ing of the gas (see Ferland 1996 for more details). 

Photoionization calculations are often parameterized in 
terms of the ionization parameter. The reason for this is the 
similarity of the results obtained for different combinations 
of Q, e, and nu giving rise to the same (U) in equation (Q) 
(e.g., Evans & Dopita 1985; McCall et al. 1985; McGaugh 
1991; Stasihska & Leitherer 1996). By analogy, it is conve- 
nient here to fix M, and hence Q(t') (eq. |g h a nd use as a 
variable the value of (G)(t') at t' = (eq. jio)). We refer 
to this quantity in the following as the 'zero-age effective 
ionization parameter', noted 



{u ) = (G)(0). (ii) 

Our results do not depend at all on the choice of M* (§3.1), 
except through the fact that any M* imposes a maximum 
(Go) at fixed hn (eqs. (jl and Jul). For all but very high 
(Go), therefore, we adopt M*= 3 x 10 4 M©, corresponding 
to Q(0) si9x lO^s -1 for mu = 100 Mq and Z = Z Q (i.e., 
the equivalent of roughly 80 07 V stars; Vacca 1994). This 
is typical of the rate of ionizing photons in giant H n regions 
in the Milky Way and other nearby galaxies (e.g., Kennicutt 
1984; for (G ) £ 10~ 15 , we adopt M,= 1 x 10 6 Mq). For 
given Sh and (G), equations (^) and ( |lo| ) allow us to com- 
pute e (we adopt olb = 2.59 x 10~ 13 cm sec -1 , appropriate 
for electronic temperatures T c ~ 10 4 K; Osterbrock 1989). 

Throughout this paper, we assume that galaxies are ion- 
ization bounded. This is suggested by the upper limit of 
3 per cent on the fraction of the ionizing radiation escap- 
ing from nearby starburst galaxies (Leitherer et al. 1995; 
see also Giallongo, Fontana & Madau 1997). We stop the 
photoionization calculations when the electron density falls 
below 1 per cent of the hydrogen density or if the temper- 
ature falls below 100 K. A model Hn region, as defined by 
equation (^), consists typically of an inner zone containing 
highly-ionized species (e.g., N +2 , +2 ), an intermediate zone 
(N + , + ), and an outer zone containing significant H° and 
other neutral and low-ionization species (O , S + ; see for ex- 
ample Fig. 1 of Evans & Dopita 1985). The relative sizes of 
the different zones depend on the model parameters. 

2.4 Emergent Luminosity 

We have now specified our model for computing the lu- 
minosity L~^(t) produced at the time t by a star- forming 
galaxy as defined by equation (^). In summary, the main 
adjustable parameters describing the luminosity S\ emitted 
by the stars are: the upper cutoff mass of the IMF, mu; 
the star formation rate, ip; and the metallicity, Z. For the 
transmission function of the photoionized gas (of effec- 
tive metallicity Z = Z), the main adjustable parameters are: 
the zero-age effective ionization parameter, (Go); the effec- 
tive gas density, hn; and the effective dust-to-heavy element 
ratio, £d- The other main parameter is the age of the galaxy, 
t. In practice, we compute the transmission function T£ of 
the gas by subdividing star formation into units of mass M* . 
This amounts to dividing ip(t — t') and multiplying S\(t') by 
M* in equation (|l|) and then using equation (g) to evaluate 
Q(t'). 

Until now, we have ignored the absorption of photons 
emanating from the ionized gas (and those emitted by old 
stars) by dust in the 'neutral ISM' before they escape from 
the galaxy. A detailed modelling of this effect is not nec- 
essary if the only observables under interest are the lu- 
minosities of optical emission lines, and if at least two H- 
recombination lines are available. In this case, the observed 
line luminosities can be corrected in a straightforward way 
for absorption by dust in the neutral ISM. The reason for 
this is that the ionized gas always produces H-recombination 
lines in ratios corresponding to dust-free case B recombi- 
nation, whether or not dust is present (Hummer & Storey 
1992; Ferland 1996; Bottorff et al. 1998; see also §4 below). 
The effect of dust is to reduce the ionization rate and hence 
the luminosities of all H-recombination lines by a similar 
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amount. The departure of H-recombination line ratios (e.g., 
Hq/H/3) from dust-free case B recombination is therefore an 
efficient probe of dust in the neutral ISM of a galaxy. In ad- 
dition, since most known extinction curves (e.g., Milky Way, 
LMC, SMC, starburst galaxies) are similar at optical wave- 
lengths, the absorption of other optical emission lines (e.g., 
[On], [Oiii], [Nil], [Sn]) can be inferred reliably from that 
of H-recombination lines. In §3.1, we use this remarkable 
property to calibrate our model for Lf- 

A more refined description of the absorption by dust 
in the neutral ISM is required to compute simultaneously 
the line and continuum emission from galaxies. In nearby 
starburst galaxies, the attenuation inferred from the Hq/H/3 
ratio is typically higher than that inferred from the spectral 
continuum (e.g., Fanelli, O'Connell & Thuan 1988; Calzetti, 
Kinney & Storchi-Bergmann 1994; Calzetti 1997). Recently, 
Chariot & Fall (2000) have developed a simple model for the 
absorption of starlight by dust in galaxies, based on an ide- 
alized description of the main features of the ISM. Chariot 
& Fall show that the finite lifetimes of the dense clouds in 
which stars form are a key ingredient for resolving the ap- 
parent discrepancy between the attenuation of line and con- 
tinuum photons in nearby starburst galaxies. Their model 
accounts for all the available observations of these galaxies, 
including the ratio of far-infrared to ultraviolet luminosities, 
the Ha/H/3 ratio, the Ha equivalent width, and the ultravi- 
olet spectral slope. We use this prescription here to compute 
the absorption by dust in the neutral ISM. 

The luminosity per unit wavelength L\{i) emerging at 
time t from a star-forming galaxy can be inferred from the 
expression of L~^(t) in equation (Fy). We write 



L x (t)= / dt'^{t-t')S x {t')T+{t')Tl{t'), 



(12) 



where T x (t') is the transmission function of the neutral ISM. 
Chariot & Fall (2000) show that the observed relationships 
between the various integrated spectral properties of nearby 
starburst galaxies are well reproduced by the simple recipejj] 

T A °(t') = exp [-f A (t')] , (13) 



TV 



(A/5500 A) 
firv (A/5500 A) 



-0.7 

-0.7 



for t' < 10 7 yr, 
for H > 10 7 yr . 



(14) 



The wavelength dependence of the 'effective absorption' 
curve t\ is tightly constrained by an observed relation be- 
tween ratio of far-infrared to ultraviolet luminosities and 
ultraviolet spectral slope. The age 10 7 yr in the above ex- 
pression, which coincides with the value of t OB in our model 
(§2.2), corresponds to the typical timescale for young stars 
to disrupt their dense 'birth clouds' or migrate away from 
them into the 'ambient ISM'. Here, we have introduced for 
convenience the fraction \x of the total effective absorption 
optical depth in the neutral ISM (birth clouds plus ambi- 
ent ISM) contributed by the ambient ISM. Chariot & Fall 
(2000) show that /j, « 1/3 approximates closely the observed 
relation between Ha/H/3 ratio and ultraviolet spectral slope 

t The function denoted here by T?(t') is equivalent to the func- 
tion T x (t') defined by equation (6) of Chariot & Fall (2000) for 
/ = 0.0 in their notations. 



for starburst galaxies. To account for the scatter about this 
relation, we consider below values ranging from fi — 1 (i.e., 
all the radiation from the H II regions leaks into the ambient 
ISM) to /i = 1/5. We note that, for fv = 0, equations ([!]) 
and ( p^[ ) are equivalent, and L\{t) = L^(t). 

We compute the luminosity L\(t) as the sum 



L\(t) = L c x {t) + S(X — Xi)Li(t) , 



(15) 



where L\(i) is the luminosity per unit wavelength of the 
continuum radiation emerging from the galaxy (including 
the gas recombination continuum), Li(t) is the luminosity 
of the line emerging at the wavelength Ai, and 8 is the Dirac 
delta function. The equivalent width of an emergent line, 
noted Wi(t), is the difference between the emission equiva- 
lent width [i.e., the ratio of Li(i) to the continuum luminos- 
ity near Ai] and the stellar absorption equivalent width. In 
star-forming galaxies, stellar absorption is especially impor- 
tant for H-recombination lines (see §2.2). Finally, we relate 
the emergent luminosity of an emission line to the star for- 
mation rate by the 'line efficiency factor' 



Vi (t) = Li(t)/ij(t). 



(16) 



In the following, we are primarily interested in the depen- 
dence of r/Ha , V[o ii] i an d V[o in] on the star-formation param- 
eters of galaxies. 



3 COMPARISON WITH OBSERVATIONS 

In this section, we use observations of nearby star-forming 
galaxies to calibrate our model. Our goal is to identify 
models that account for the typical properties, scatter, and 
trends seen in the observations. In §3.1 below, we first de- 
scribe the specific influence of each parameter in the model 
on the ratios of prominent optical emission lines. This allows 
us to identify how the Ha, [On], and [O ill] efficiency fac- 
tors are related to observed line ratios. This analysis, which 
involves line luminosities only, does not depend on our as- 
sumptions about the absorption by dust in the neutral ISM 
(§2.4). We must check, however, that the model also repro- 
duces whole (line plus continuum) spectral energy distribu- 
tions of nearby star-forming galaxies. We show this in §3.2, 
where we then have to specify the absorption by dust in the 
neutral ISM. 



3.1 Emission-Line Luminosities 

To calibrate our model, we compare it with observations 
of a sample of nearby star-forming galaxies, for which the 
fluxes of the Ha, H/3, [On]A3727, [Om]A5007, [Nn]A6583, 
and [Sll]AA6717, 6731 emission lines are available. The sam- 
ple, which we have compiled from the literature, includes 92 
non-Seyfert galaxies spanning a wide range of morphological 
types, from Sab to Irr, and a wide range of absolute B-band 
magnitudes, — 16 ^ Mb —22 (we adopt a Hubble constant 
Ho = 70 km s - 1 Mpc~ 1 ) . About a third of these (32 galaxies) 
are normal spiral and irregular galaxies from the sample of 
Kennicutt (1992b). Another third (33 galaxies) are starburst 
galaxies selected by Calzetti, Kinney & Storchi-Bergmann 
(1994) from the International Ultraviolet Explorer (IUE) at- 
las of Kinney et al. (1993), for which optical-line fluxes are 
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available from Storchi-Bergmann, Kinney & Challis (1995) 
and McQuade, Calzetti & Kinney (1995). The remaining 
third (27 galaxies) are low-metallicity, 'H II galaxies' (in- 
cluding isolated extragalactic H II regions and blue compact 
dwarf galaxies) from the sample of Stasiriska & Leitherer 
(1996). These include 17 galaxies with secure Ha measure- 
ments from the Terlevich et al. (1991) catalogue and 10 
galaxies observed by Izotov, Thuan & Lipovetzky (1994). 
Most of the normal spiral and starburst galaxies in our sam- 
ple have Ha equivalent widths in the range 10-150 A. The 
low-metallicity, H II galaxies have typically much higher Ha 
equivalent widths, in the range 300-800 A. 

We correct the flux ratios of all lines for absorption 
by dust in the neutral ISM on the basis of the observed 
Ha/H/3 ratios (§2.4). This requires us to first correct the ob- 
served Ha and H/3 fluxes for stellar absorption. Calzetti et 
al. (1994) and Izotov et al. (1994) derive H/3 stellar absorp- 
tion equivalent widths for all the galaxies in their samples 
using the procedure outlined by McCall et al. (1985), that is 
based on a simultaneous fitting of the observed Ha, H/3, and 
H7 fluxes. Also, Kennicutt (1992b) shows that a mean H/3 
absorption equivalent width of —5 A is appropriate for the 
galaxies in his sample, and Masegosa et al. (1994) suggest a 
mean H/3 absorption equivalent width of —2 A for the low- 
metallicity galaxies of the Terlevich et al. (1991) catalogue 
(this is also roughly the mean value for the galaxies of the 
Izotov et al. 1994 sample). To correct the observed Ha/H/3 
ratio of each galaxy for stellar absorption, we adopt the same 
absorption equivalent width for Ha as for H/3 (see §2.2 for a 
justification). We then evaluate the effective absorption by 
dust in the neutral ISM by comparing the resulting Ha/H/3 
ratio to the dust-free case B recombination value of 2.86 
(appropriate for electronic densities n c 10 4 cm~ 3 and tem- 
peratures T c w 10 4 K; Osterbrock 1989). Finally, we dered- 
den the relative ratios of the [On], H/3, [Om], Ha, [Nil], 
and [S 11] lines using the Calzetti (1997) effective absorption 
curve (which is similar to the Milky Way, LMC, and SMC ex- 
tinction curves at optical wavelengths). Neglecting possible 
anisotropies, we equate line flux ratios to the correspond- 
ing luminosity ratios. Figs. 1, 2, and 3 show £jo m ] 
as a function of L+ n] /L+ m] , L^ Sn] /L+ a , and L+ u] /L+ 
respectively, for this sample. Evidently, the galaxies follow 
well-defined relations in these diagrams. The typical mea- 
surement errors indicate that the scatter about the mean 
relations is real. 

We now use these observations to constrain the param- 
eters in our model. Our goal is to identify models that ac- 
count for the typical properties, scatter, and trends seen in 
the sample in order to relate the Ha, [On], and [O in] effi- 
ciency factors to observed line ratios. The influence of each 
parameter on observable quantities can best be explored by 
keeping all other parameters fixed at 'standard' values. For 
simplicity, we take for the moment the star formation rate 
ip to be constant. In this case, the age t of the stellar pop- 
ulation should be regarded as the effective age of the most 
recent burst of star formation. After some experimentation, 
we adopted the following standard parameters: 
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Figure 1. [O in] /H/3 ratio plotted against [Oll]/[Om] ratio. The 
absorption by dust in the neutral ISM and the stellar H/3 absorp- 
tion are not included. The data points are from the samples of 
galaxies and individual Hn regions discussed in §3.1 and are re- 
peated in all panels (squares: normal spiral and irregular galaxies; 
stars: starburst galaxies; crosses: Hn galaxies; and dots: Hn re- 
gions; typical measurement errors for galaxies are indicated in the 
upper left panel). In each case, the line shows the effect of vary- 
ing one parameter from the lower end of the range (square) to the 
standard value (circle) to the upper end of the range (triangle), 
with all other parameters fixed at their standard values (eq. [1171] ) : 
(a) effective starburst age, t = 1 X 10 6 , 3 X 10 s , and 3 X 10 9 yr](b) 
upper cutoff mass of the IMF, m v = 40, 100, and 130 M Q ; (c) 
zero-age effective ionization parameter: log((f7o)) = —3.0, —2.5, 
and -2.0; (d) effective metallicity: Z = 0.2Z Q , Zq, and 2Z Q ; (e) 
effective dust-to-heavy element ratio: £d = 0.1, 0.3, and 0.5 (f) 
effective gas density: hu = 3, 30, and 300 cm -3 . 



UK 



0.3 

30 cm" 



(17) 



While these values are not the result of a rigorous opti- 
mization procedure, they do enable the standard model to 
match roughly the observed typical (i.e., median) properties 
of the galaxy sample. We note that £d ~ 0.3 is of the order 
of the dust-to- heavy element ratio expected in Hn regions 
(§2.3). Also, fiu ~ 30cm~ 3 is typical of the gas density in 
Galactic and extragalactic H 11 regions (10-100 cm -3 ; Oster- 
brock 1989). For this hn, log((J7o)) ~ -2.5 implies e » 0.04 
(§2.4), which is typical of the volume-filling factor in giant 
Hn regions in the Milky Way and other nearby galaxies 
(e.g., Kennicutt 1984). 

We now compare the predictions of our model with the 
observations described above. Each panel in Figs. 1-3 shows 
the effect of increasing and decreasing one parameter with 
respect to its standard value with the others held fixed. In 
Fig. 4, we show the influence of these variations on the Ha, 
[0 11], and [O in] efficiency factors r/ua , Wo 11] , an d V[o ml com- 
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Figure 2. [O in] /H/3 ratio plotted against [Sn]/Ho ratio. The 
absorption by dust in the neutral ISM and the stellar Ha and H/3 
absorptions are not included. The data points are from the sam- 
ples of galaxies and individual Hn regions discussed in §3.1 and 
are repeated in all panels (squares: normal spiral and irregular 
galaxies; stars: starburst galaxies; crosses: Hn galaxies; and dots: 
H II regions; typical measurement errors for galaxies are indicated 
in the upper left panel). In each case, the line shows the effect of 
varying one parameter from the lower end of the range (square) 
to the standard value (circle) to the upper end of the range (tri- 
angle), with all other parameters fixed at their standard values 
(cq. 117 ). The models are the same as in Fig. 1. 
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Figure 3. [Om]/H/3 ratio plotted against [Nii]/[Sii] ratio. The 
absorption by dust in the neutral ISM and the stellar H/3 ab- 
sorption are not included. The data points are from the samples 
of galaxies and individual Hll regions discussed in §3.1 and are 
repeated in all panels (squares: normal spiral and irregular galax- 
ies; stars: starburst galaxies; crosses: H II galaxies; and dots: H II 
regions; typical measurement errors for galaxies are indicated in 
the upper left panel). In each case, the line shows the effect of 
varying one parameter from the lower end of the range (square) 
to the standard value (circle) to the upper end of the range (tri- 
angle), with all other parameters fixed at their standard values 
(cq. |17|). The models are the same as in Fig. 1. 



puted using equation ( [iq ) with f v = (i.e., for Li = Lf). 
We can summarize the role of each parameter as follows. 

Effective starburst age. Since we have assumed tp — const, 
increasing t makes rjn a , tj^q n | , and »?[oiii] larger initially, as 
massive stars with lifetimes 3-5xl0 6 yr accumulate on the 
main sequence (Fig. 4a). Stars with ages between 5xl0 6 yr 
and t OB = 1 x 10 7 yr produce less ionizing photons, while 
older stars do not contribute to the ionizing radiation in our 
model (§2). After 1 x 10 7 yr, therefore, a steady population 
of H II regions with effective ionization parameters ranging 
from (U ) to (U){t OB ) is established (eq. @). The rising 
contribution of gas with (U)(t) < (Uo) at early ages has a 
small but significant effect on the line ratios in Figs. la-3a, 
which can be understood from our discussion of (Uo) below. 

Upper cutoff mass of the IMF. Increasing mu makes rjHa, 
rjyo „] , and ??[o m ] larger because, near the cutoff, more mas- 
sive stars produce more ionizing photons. This effect be- 
comes negligible for mu > 100 Mq (Fig. 46). Increasing mu 
and hence Q(t') (eq. ||) also makes the effective ionization 
parameter larger in the inner (0 +2 ) parts of the Hll regions 
(eq. |J). At fixed {Uo), however, the effective ionization pa- 
rameter in the outer, low-ionization (S^) boundaries near 
the Stromgren radius is also fixed (eq. |l0[|). Thus, increas- 
ing mu makes the [O III] /H/3 ratio larger and the [O n]/[0 m] 



ratio smaller but affects weakly the [Sll]/Htv and [Nll]/[Sn] 
ratios (Figs. 16-36; see §2.3). 

Zero-age effective ionization parameter. Increasing (Uo), 
which corresponds here to increasing e at fixed tt-h and 
Q(0), causes more gas to be concentrated in the inner, high- 
ionization parts of the Hll regions; thus the [Olll]/H/3 ratio 
and r7[oiii] increase, and the [On]/[Oni] ratio and t?[oh] de- 
crease (Figs, lc and 4c). Increasing (Uo) also increases the ef- 
fective ionization parameter near the Stromgren radius (Rs) 
and weakens the low-ionization [Sn] line (Figs. 2c and 3c). 
Since the H-column density scales roughly as Inn Rs oc (U), 
increasing (Uo) at fixed dust-to-gas ratio £dZ causes more 
ionizing photons to be absorbed by dust and hence reduces 
r/Hc (Fig. 4c). 

Effective Metallicity . Increasing Z makes the cooling by 
heavy elements more efficient and reduces the electronic 
temperature. The cooling becomes dominated by infrared 
fine-structure transitions with low excitation energies, and 
fewer electrons are capable of exciting optical lines (e.g., 
Spitzer 1978). Thus, the [O iii]/H/3 ratio and Tyron] an d 7?[om] 
decline (Figs. Id and 4d). The [Oii]/[Oiii] ratio increases 
in Fig. Id because of the efficient fine-structure cooling by 
doubly- ionized species in the inner parts of the Hll regions 
(Stasihska 1980). The [Sll]/Ha ratio, which increases ini- 
tially due to the higher S + abundance, also declines when 
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Figure 4. Ha, [On], and [Oin] efficiency factors (as defined by 
cq. normalized to the standard model of eq. Jit] ]. The ab- 

sorption by dust in the neutral ISM and the stellar Ha absorption 
are not included. In each panel, the line shows the effect of vary- 
ing one parameter from the lower end of the range (square) to the 
standard value (circle) to the upper end of the range (triangle), 
with all other parameters fixed at their standard values (eq. [11 71]). 
The models are the same as in Fig. 1. 

cooling becomes important (Fig. 2d). In Fig. 3d, the increase 
of the [Nii]/[Sii] ratio with Z follows from our inclusion of 
nitrogen from secondary processing (§2.3). We note that the 
rate of ionizing photons produced by massive stars tends to 
decrease with increasing metallicity (e.g., Fig. 4 of Garcia- 
Vargas et al. 1995). This effect accounts for roughly two 
thirds of the drop in rjua as Z increases in Fig. Ad, one third 
arising from the larger fraction of ionizing photons absorbed 
by dust as ^Z increases. 

Effective dust-to-heavy element ratio . Increasing £d de- 
pletes important coolants such as oxygen from the gas phase. 
As the cooling through the infrared fine-structure transitions 
of these elements declines, the electronic temperature in- 
creases and makes cooling through optical transitions more 
efficient (Shields & Kennicutt 1995). The [Oill]/H/3 ratio 
and Won] an d r?[om] increase only slightly in Figs, le and 
4e because the rise in T e is compensated by the depletion of 
oxygen from the gas phase. However, since S and N are not 
refractory elements, the [Sll]/Ha ratio increases substan- 
tially with £a at roughly constant [Nn]/[Sn] ratio (Figs. 2e 
and 3e). Increasing £d also causes more ionizing photons to 
be absorbed by dust and hence reduces r/Ha (Fig- 4e). 

Effective gas density. Increasing hn makes collisional de- 
excitation compete with radiative cooling. Since the infrared 
fine-structure transitions have lower critical densities for col- 
lisional de-excitation than the optical transitions, their cool- 
ing efficiency declines when increasing Sh (Stasinska 1990). 
The effect is negligible at low metallicities but more signifi- 
cant at high metallicities, where the fine-structure lines are 



the dominant coolants (Oey & Kennicutt 1993). The result- 
ing rise in T e makes the cooling through optical transitions 
more important ; hence ^ron] and ^7[o m] 

increase slightly 

(Figs. 4f). Changes in Ah have a negligible effect on r)n a 
and a small effect on the line ratios in Figs. 1/-3/. 

Other parameters . We have also computed models with 
different IMF slopes (not shown). Changes in x have a neg- 
ligible effect on the line ratios, which are determined by the 
properties of the most massive stars. For a power-law IMF, 
however, increasing x causes more mass to be locked in low- 
mass stars and reduces the line efficiency factors. We also 
tested the influence of changing the effective mass of the 
ionizing star clusters and found that models with different 
M* and fixed (Uo) are virtually identical. 

Since each parameter in our model has a specific influ- 
ence on the ratios of prominent optical emission lines in star- 
forming galaxies, we can comment on the physical origin of 
the relations defined by these observed quantities. Figs. 1-3 
show that the observed mean trends can be most naturally 
reproduced by a sequence in the effective gas metallicity, 
although a spread in the effective ionization parameter is 
also indicated. In fact, at one end of the observed relations, 
the very low [Oii]/[Oiii] ratios of the Hn galaxies with the 
smallest [Nii]/[Sii] ratios appear to require high effective 
ionization parameters in addition to low metallicities (Izotov 
et al. 1994 estimate Z ps O.IZq from a detailed spectral anal- 
ysis of these galaxies; this is close to the value Z = Q.2Zq 
of the lowest-metallicity model in Figs ld-3d). At the other 
end, the [Om]/H/3, [O ll]/[0 in], and [Nii]/[Sii] ratios of the 
most extreme normal spiral and starburst galaxies appear 
to require low effective ionization parameters in addition to 
high metallicities. As Figs. 1-3 indicate, the intrinsic scatter 
about the observed relations can arise from variations in the 
effective dust-to-heavy element ratio, the effective gas den- 
sity, and the effective age of the starburst (and possibly the 
IMF). 

The primary dependence of optical emission-line ratios 
on the metallicity and ionization parameter of the gas is 
a well-known property of individual Hn regions (e.g., Bald- 
win, Phillips & Terlevich 1981; Evans & Dopita 1985; McCall 
et al. 1985; Campbell 1988; McGaugh 1991; Shields & Ken- 
nicutt 1995; Bresolin et al. 1999). To illustrate the parallel 
between galaxies and individual H II regions, we also show in 
Figs. 1-3 observations of 68 H II regions in nearby spiral and 
irregular galaxies from McCall et al. (1985). The observed 
relations for Hn regions are similar to those for galaxies. 
This similarity, which justifies a posteriori our adoption of 
effective parameters to describe the nebular properties of 
whole galaxies, arises presumably from the limited range of 
metallicities and ionization parameters of the gas in most 
star- forming galaxies (e.g., Zaritsky, Kennicutt & Huchra 
1994; Kobulnicky et al. 1999). The rough coincidence be- 
tween starburst and normal spiral galaxies in Figs. 1-3 fur- 
ther suggests that stochastic bursts of star formation do not 
strongly influence the global gas parameters of galaxies. 

Fig. 2 does reveal a difference between the nebular prop- 
erties of galaxies and those of individual H II regions. At fixed 
[O iii]/H/3 ratio, galaxies tend to have higher [S ll]/H« ratios 
than individual Hn regions. This enhancement in the rela- 
tive luminosity of the low-ionization [S n] line is thought to 
be a signature of the diffuse ionized gas in galaxies (e.g., 
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Lehnert & Heckman 1994; Martin 1997; Wang et al. 1997). 
It is usually associated with enhanced [Ol]A6300, [Nil], and 
[Oil] emission and reduced [Oiii] emission. Fig. 2c con- 
firms that the enhanced [Sll]/Ha ratios of galaxies rela- 
tive to individual Hn regions can be interpreted in terms 
of a lower effective ionization parameter at fixed effective 
metallicity. We note, however, that the effective dust-to- 
heavy element ratio in the ionized gas has also a major 
influence on the [Sll]/Ha ratio (Fig. 2e). While decreasing 
(Uo) makes the [Oiii]/H/3 ratio smaller, increasing £d allows 
models to reproduce the observations of galaxies with both 
large [S ll]/Ha and [O mJ/H/3 ratios (the effects of £d on the 
[Nil] /Ha and [Ol]/Ha ratios are similar to those on the 
[Sii]/Hq and [Oiii]/H/3 ratios, respectively). It is worth re- 
calling that, in our model, the details of the Hn regions and 
the diffuse gas ionized by each stellar generation in a galaxy 
are subsumed in our description of the gas by effective pa- 
rameters. 

One of the most important novelties of our model is 
that it allows us to relate line efficiency factors to observed 
optical-line ratios. This has fundamental implications for 
SFR estimates, as can be illustrated by the following ex- 
ample. Figs. 1 and 4 show that increasing (Uo) or decreas- 
ing Z relative to our standard model have similar effects on 
the observed [O n]/[0 ill] and [O m]/H/8 ratios, but opposite 
effects on J?[on]- Thus, for fixed [Oii]/[Oiii] and [Om]/H/3 
ratios, there can be a spread of at least a factor of several 
in J7[o n] , depending on the effective metallicity and effective 
ionization parameter of the gas. Fig. 3 shows, however, that 
the [Nii]/[Sii] ratio can help us distinguish between varia- 
tions in (Uo) and Z. Therefore, by combining information 
about the [OlIl]/H/3, [O ll]/[0 ill], and [Nii]/[Sii] ratios, we 
can reduce the uncertainties in Tjron] an d hence better con- 
strain the star formation rate. We exploit this property of 
our model in §4 below, where we construct estimators of 
the star formation rate for different assumptions about the 
available spectral information. 

3.2 Spectral Energy Distributions 

Up to now, we have compared our model with observations 
only in terms of emission-line luminosities. Since we com- 
pute in a consistent way the emission from stars and gas, 
it is important to check that the model can also reproduce 
whole (line plus continuum) spectral energy distributions of 
observed galaxies. This requires us to consider two complica- 
tions: the contamination of H-recombination lines by stellar 
absorption; and the absorption of line and continuum pho- 
tons by dust in the neutral ISM. In §2, we described how 
these effects are accounted for in our model. 

We now compare our model with observed spectral en- 
ergy distributions of various types of nearby star-forming 
galaxies. The integrated spectra at wavelengths 3650-7100 A 
of the galaxies in the Kennicutt (1992b) sample used in §3.1 
are available from Kennicutt (1992a). To avoid using the 
spectrum of any arbitrary galaxy as reference, we extract 
observed 'templates' from this sample by averaging the spec- 
tra of a few galaxies of the same morphological type with 
roughly similar emission-line properties. Figs. 5, 6, and 7 
show such spectral templates for Sb, Sc, and (Magellanic- 
irregular) Sm/Im galaxies, respectively. Each spectrum is 
an average over three galaxies with comparable [Oii]/[Oiii] 
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Figure 5. Spectral fit of an observed Sb galaxy template. The 
template (thin line) is an average over three Sb galaxies with com- 
parable [On] /[Oiii] and Ho/H/3 ratios in the Kennicutt (1992a) 
atlas (NGC 1832, NGC 3147, and NGC 3627). The model (thick 
line) has a exponentially declining star formation rate (timescale 
5 Gyr) and an age t ~ 9 Gyr. The gas parameters are: log(([/o)) = 
—3.0, Z = 1.5Zq, £ d = 0.3, and fy = 1.0. The upper four panels 
show the most prominent lines on a larger wavelength scale (a 
mark indicates the scale A A = 30 A). 



and Ha/H/3 ratios. For reference, the [O n]/[0 ill] ratios and 
the Ha and H/3 equivalent widths of the Sb, Sc, and Sm/Im 
templates are L[on]/£[oin] ~ 3.3, 2.6, and 1.4, Whu ~ 14, 
26, and 82 A, and W H /3 ~ -0.7, 3.0, and 14 A, respectively. 
The spectral resolution is about 800. 

To reproduce these observed templates, we use various 
spectral features to adjust the different parameters in our 
model. The equivalent widths of emission lines and the stel- 
lar continuum (in particular, the 4000 A discontinuity) con- 
strain the ratio of young (blue) to old (red) stars in the 
model and hence the star formation history (e.g., Kenni- 
cutt 1983; Bruzual & Chariot 1993). In principle, these con- 
straints would depend on the unknown metallicity and the 
unknown dust content of the galaxies. However, since the 
star formation history fixes the stellar Ha and H/3 absorp- 
tions, the Ha/H/3 ratio provides an independent measure 
of the absorption by dust. Then, using the [Nii]/[Sii] and 
[Oii]/[Oiii] ratios, we can adjust the effective metallicity 
and the effective ionization parameter of the gas (Figs, lc 
and 3d). Finally, for fixed values of these parameters, the 
[S ll]/Ha ratio constrains the dust content of the ionized gas 
(Fig. 2e). Therefore, we can tune in a methodic way the 
parameters of our model to reproduce the observations. 

To select models that are consistent with the observa- 
tions, we explore only a coarse grid of parameters (see Ta- 
ble |l] below). Attempting to fine tune the parameters might 
be of limited validity because the observed spectral tem- 
plates were built from individual galaxies that may have 
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Figure 6. Spectral fit of an observed Sc galaxy template. The 
template (thin line) is an average over three Sc galaxies with com- 
parable [Oii]/[Oiii] and Hoj/H/3 ratios in the Kennicutt (1992a) 
atlas (NGC 2276, NGC 2903, and NGC 6181). The model (thick 
line) has a exponentially declining star formation rate (timcscale 
15 Gyr) and an age t sa 8 Gyr. The gas parameters are: 
log«£/ )) = -3-0, Z = Z Q , £ d = 0.1, and r v = 0.8. The upper 
four panels show the most prominent lines on a larger wavelength 
scale (a mark indicates the scale A A = 30 A). 
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Figure 7. Spectral fit of an observed Sm/Im galaxy template. 
The template (thin line) is an average over three (Magellanic- 
irregular) Sm/Im galaxies with comparable [Oll]/[Om] and 
Ha/H/3 ratios in the Kennicutt (1992a) atlas (NGC 4449, 
NGC 4485, and NGC 4670). The model (thick line) has a constant 
star formation rate and an age t sa 2.5 Gyr. The gas parameters 
are: log«£?o» = -3.0, Z = O.2Z , £ d = 0.1, and f v = 0.1. 
The upper four panels show the most prominent lines on a larger 
wavelength scale (a mark indicates the scale A A = 30 A). 



slightly different properties. We take the star formation rate 
ip(t) to be either constant or exponentially declining. In 
practice, the results are mostly sensitive to the ratio of the 
current SFR to the average past SFR. We also explore varia- 
tions in the zero-age effective ionization parameter, (t/o), the 
effective metallicity, Z, the effective dust-to-heavy element 
ratio in the ionized gas, £d, and the effective absorption op- 
tical depth at 5500 A in the neutral ISM, fy (eqs. jl3| and 
Jl4|l). We fix the other parameters at their standard values, 
i.e., p. = 1/3 (eq. (m)), m v = 100 M Q , and n H = 30cm -3 
(eq. Figs. 5, 6, and 7 show models selected in this 

way that reproduce the observed Sb, Sc, and Sm/Im spec- 
tral templates. In each case, the continuum spectral energy 
distribution and the luminosities of all the optical emission 
lines in the models agree well with the observations. Given 
the idealizations of the models, we consider these results 
as very satisfactory. The line strengths in the models are 
also quantitatively close to those quoted above for the ob- 
served spectral templates. For reference, the model Sb, Sc, 
and Sm/Im spectra, respectively, have L[on]/£[Om] ~ 2.9, 
2.5, and 1.6, Wn a « 14, 30, and 89 A, and W Kf3 ~ -0.4, 2.3, 
and 14 A. The contribution to these equivalent widths by 
stellar absorption is —5.2, —6.2, and —7.3 A, respectively. 

These spectral fits constrain simultaneously the stars, 
gas, and dust parameters of the galaxies. We find that the 
ratio of the current SFR to the average past SFR favored by 
the observations ranges from 0.4 for the Sb galaxy template 
to 0.8 for the Sc galaxy template to 1.0 for the Sm/Im galaxy 
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Figure 8. Spectral energy distributions of the same models as in 
Figs. 5—7 shown with extensions at ultraviolet and near-infrared 
wavelengths. 
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template. This is consistent with the expected dependence of 
this ratio on morphological type (e.g., Kennicutt, Tamblyn 
& Congdon 1994). Furthermore, the emission- line ratios of 
the Sb and Sc galaxy templates are well reproduced by mod- 
els with roughly solar effective metallicity and fy » 1, while 
the (Magellanic-irregular) Sm/Im galaxy template is better 
fitted by a model with both low effective metallicity and 
dust content. All three spectra favor the same zero-age ef- 
fective ionization parameter, (Uo) ~ —3.0, but slightly dif- 
ferent effective dust-to-heavy element ratios in the ionized 
gas. Of course, these are only rough constraints on the star 
formation properties of the observed galaxies, as we have 
deliberately avoided fine tuning the model parameters (the 
uncertainties in the derivation of the physical parameters 
of galaxies will be addressed in §4). In addition, the ob- 
served templates in Fig. 5-7 are representative of only some 
Sb, Sc, and Sm/Im galaxies (see above). Thus, these results 
should be taken merely as illustrative of the ability with 
our model to constrain simultaneously the stars, gas, and 
dust parameters of galaxies using integrated spectral energy 
distributions. 

We have described here the optical properties of our 
model for the emission from star-forming galaxies. In reality, 
we compute spectral energy distributions of the line and con- 
tinuum emission over the entire range of wavelengths from 
91 A to 5 fim (redward of this limit, the emission from poly- 
cyclic aromatic hydrocarbons, which are not included in the 
model, becomes important). As an example, Fig. 8 shows 
extensions of the model spectra in Figs. 5-7 at ultraviolet 
and near-infrared wavelengths. The prominent emission lines 
outside the optical window (e.g., Mgn at A = 2798 A and 
[S in], [Fen], Pa at near-infrared wavelengths) offer comple- 
mentary diagnostics of the star formation parameters, which 
we plan to study in future work. Finally, it is worth recalling 
that the prescription we have adopted to describe the ab- 
sorption by dust in the neutral ISM ensures that our model is 
also consistent with the observed ultraviolet and far-infrared 
properties of nearby starburst galaxies (see §2.4 above). 



4 OBSERVATIONAL CONSTRAINTS ON THE 
STARS, GAS, AND DUST PARAMETERS 
OF GALAXIES 

We now have a model for relating the integrated spectral 
properties of galaxies to physical parameters such as the star 
formation rate, the metallicity, and the effective absorption 
optical depth of the dust. In practice, whole spectral energy 
distributions of the type analyzed in §3.2 are not always 
available to constrain these parameters in galaxies. However, 
the tight correlations between the different integrated spec- 
tral properties of nearby star-forming galaxies suggest that 
even partial spectral information may be useful to constrain 
physical parameters (Figs. 1-3 above; see also Figs. 2 and 
3 of Chariot & Fall 2000). Here we analyze these observed 
relations with our model to construct estimators of the star 
formation rate, the gas-phase oxygen abundance, and the 
effective absorption optical depth of the dust in galaxies for 
various assumptions about the available spectral informa- 
tion. 

Our first step is to select a comprehensive set of mod- 
els that account for the full range of observations of nearby 



star-forming galaxies. We use the observed emission-line ra- 
tios of the galaxies in Figs. 1-3 to identify the appropriate 
ranges in zero- age effective ionization parameter, {Uo}, effec- 
tive metallicity, Z, and effective dust-to-heavy element ratio, 
£d- Since the effective gas density has only a weak influence 
on line luminosities (Figs. 1/-4/), we fix this parameter at its 
standard value of equation (|rij), i.e., hn — 30 cm -3 . We also 
consider variations in the total effective absorption optical 
depth of the dust in the neutral ISM, fv, and the fraction of 
this contributed by the ambient ISM, [i (eq. These are 

constrained by the observed relations between ratio of far- 
infrared to ultraviolet luminosities, Ha/H/3 ratio, and ultra- 
violet spectral slope in nearby starburst galaxies (Chariot 
& Fall 2000). Finally, we consider different star formation 
histories (as before, we assume that stars and gas have the 
same metallicity, i.e., Z — Z). We take the star formation 
rate ip(t) to be either constant or exponentially declining. In 
each case, we compute models at ages ranging from 10 7 yr 
to 10 10 yr. For simplicity, we adopt in all models a Salpeter 
IMF truncated at 0.1 and 100 Mq (the line efficiency fac- 
tors in Table ^ below would be a factor of about 3.4 times 
smaller for a Scalo 1986 IMF with the same cutoff masses). 

In Table [l], we report the parameters of the models 
found to reproduce the observations of nearby star-forming 
galaxies. Most but not all combinations of these parameters 
are consistent with the observations. In particular, the ob- 
served emission-line ratios of nearby galaxies (Figs. 1-3) do 
not favor combinations of both either very low or very high 
Z and (Uo). In Fig. 9, we show the emission-line ratios of 
the models with the allowed combinations of Z, {Uo), and 
| d in Table |l[ As in Figs. 1-3, we do not include in Fig. 9 
the absorption by dust in the neutral ISM and the stellar 
Ha and H/3 absorptions. Since age has a negligible effect in 
such diagrams (especially for t > 10 7 yr; see Figs. la-3a), 
we plot only models with constant ip(t) at the standard age 
t = 3 x 10 8 yr (eq. 0). Fig. 9 shows that the models in 
Table |l| account well for the typical properties, scatter, and 
trends seen in the observations. The models do not reach the 
extreme [Nii]/[Sii] ratios of some Hn and starburst galax- 
ies, suggesting that our simple prescription for secondary 
nitrogen production may break down at very low and very 
high Z (§2.3; we note that the largest observed [Nii]/[Sii] 
ratios in Fig. 9c could result from a decrease in the S/O 
abundance ratio at high metallicity; see Diaz et al. 1991). 
We have checked that the predictions of our model at these 
extreme effective metallicities do not affect our conclusions 
below. The models with various fy and fi in Table |l| are 
also consistent, by construction, with the observed ratios of 
far-infrared to ultraviolet luminosities, Ha/H/3 ratios, and 
ultraviolet spectral slopes of nearby starburst galaxies in 
Figs. 2-4 of Chariot & Fall (2000; see §2.4 above). Therefore, 
the models in Table |l| enable us to interpret the observed 
spectral properties of nearby star-forming galaxies in terms 
of stars, gas, and dust parameters. 

A striking property of these models is that they span 
vast ranges of Ha and [O n] efficiency factors, r/u a = 0.05- 
2.8 x 10 15 erg g" 1 and J7 [0n ] = 0.02-3.6 x 10 15 erg g" 1 
(for Whc and Won] both larger than 5 A). Only parts 
of these ranges arise from variations in fy and [i in Ta- 
ble |l[ For fv = 0, changes in the ionized-gas parameters 
alone {Z, {Uo), and fa) account for factors of 3.5 and 14 
spreads in ijh q and r?[on], respectively. The range of stel- 
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Table 1. Parameters of the models reproducing the observed 
integrated spectral properties of nearby star-forming galaxies. 

log(({/ >) -3.0, -2.5, -2.0, -1.5 
Z/Zq 0.2, 0.5, 0.75, 1.0, 1.5, 2.0 

| d 0.1, 0.3, 0.5 

fy 0.01, 0.3, 0.5, 1.0, 1.5, 2.0 

li- 1 1, 3, 5 

ip(t) constant, exponentially declining 

t/yr 10 7 -10 10 (112 uneven steps) 

Two timescales are adopted for the exponentially declining star 
formation rate ip(t): 0.1 and 6.0 Gyr. All models have njj = 
30 cm -3 and a Salpeter IMF truncated at 0.1 and 100 M . Most 
but not all combinations of the above parameters are allowed by 
the observations (see text for detail). 
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Figure 9. Emission-line ratios of the models in Tabic [j]. The 
absorption by dust in the neutral ISM and the stellar Ha and H/3 
absorptions are not included. Only models with constant i/j(t) are 
shown at the standard age t = 3 X 10 8 yr (eq. Jl?] ; crosses). The 
data points in (a), (ft), and (c) are the same galaxies as in Figs. 1, 
2, and 3, respectively (all types plotted as circles here). 

lar Ha absorption strengths arising from variations in the 
star formation history contributes an additional factor of 
1.8 spread in 7/Ha- These results put into perspective the ap- 
proximate efficiency factors often adopted to estimate star 
formation rates in galaxies, e.g., r]u a ~ 2.0 x 10 erg g - 
and 7][ n] ~ 1-2 X 10 15 erg g" 1 (Kennicutt 1998). Our model 
allows us to assess that the actual uncertainties in SFR es- 
timates based on either the Ha line or the [O n] line alone 
can be as high as several decades. 

Since we know how tjaa and ?7[oii] are related to ob- 
served optical-line ratios (§3.1), we can exploit this informa- 
tion to improve SFR estimates in galaxies. For the sake of 
generality, we denote here by A the quantity to be estimated 
(e.g., r/Ha, Won])- Given a set of N spectral features x k (e.g., 
line ratios, equivalent widths), we can always least-squares 
fit to the models in Table [j] an expression of A of the form 



log .4: 



(to 



N 

£ 

fe=i 



a k log x k . 



(18) 



where the a k are the fitting coefficients. It is convenient to 
introduce the logarithmic difference A m between the true 
and fitted A values of the mth model, 
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a k logx k ,r, 



(19) 



The probability distribution of the models being unknown, 
the maximum absolute value of A m over all the models in 
Table 111 provides a lower limit on the accuracy of expres- 
sion (181). By experimenting with several different combina- 



tions of spectral features x k , we can identify the one leading 
to the highest accuracy and hence the best observational 
estimate of A. 

We have used this procedure to derive estimators of not 
only the Ha and [O n] efficiency factors, but also the gas- 
phase oxygen abundance and the effective absorption optical 
depth of the dust in star-forming galaxies. In Table we 
present our results for several examples of available spectral 
information (cases A-G). The observations are assumed not 
to be corrected for dust and stellar (H-Balmer) absorption. 
In each case, we express the efficiency factor of the most 
accurate SFR tracer (r/na , ?/Ha+ [n n] , or ri[ n] ) , the gas-phase 
oxygen abundance (O/H), and the effective dust-absorption 
factor at 5500 A (e~ Tv ) as functions of the available spectral 
features. For reference, we also give the expression of the line 
efficiency factor for fv = (noted r]^ a , ??hc«+[Nii]; or V°on])- 

The accuracy of an estimator in Table depends on its 
regime of application. In Fig. 10, we plot for each estimator 
the range in A m , as defined by equation ([jl]), as a func- 
tion of the estimated quantity. The potential (logarithmic) 
error in an estimate of r\i, O/H, or e~ TV - ignoring observa- 
tional errors - is the maximum absolute value of A m at the 
abscissa corresponding to that estimate (we note that, for 
e~ TV ~ 1, even small values of |A m | can translate into large 
errors in fy). The maximum absolute value of A m over the 
full range in the estimated quantity is the maximum error 
associated to the estimator. To avoid the undue influence 
of large tails arising from a few discrepant models in the 
|A m | distributions, we define the 'uncertainty' pertaining to 
each estimator in Fig. 10 as the 99th percentile range in 
|A m | over the full range in the estimated quantity. We re- 
port this uncertainty in Table ^ and use it below to quantify 
the accuracy of an estimator (for reference, we also indicate 
in parentheses in Table ^ the full, 100th percentile range in 
|A m |). 

The results of Table H demonstrate the merit of spec- 
tral information for estimating the star formation rate, the 
metallicity, and the absorption by dust in a galaxy. In the 
best-constrained case where the [On], H/3, [Om], Ha, [Nil], 
and [Sn] lines are all observed (case A), we find that rjn a , 



O/H, and 



can be estimated to within less than a fac- 



tor of 2 (i.e., with absolute logarithmic uncertainties less 
than 0.3). This is a reduction by more than a factor of 25 
of the uncertainty in i)h relative to the case where the 
Ha line alone is used to estimate the SFR (see above). 
The reason for such an improvement is that we have in- 
cluded in the estimator the Ha/H/3, [O n]/[0 in], [O m] /H/3, 
and [Sll]/Ha ratios and the H/3 equivalent width, which 
all carry important information about the parameters that 
influence r/ua (see §3). Remarkably, a similar accuracy in 
r/ua can be achieved even if the Ha line is blended with 
the adjacent [Nil] lines (case B). In this case, however, the 
loss of the [Nll]/[Sn] ratio makes O/H estimates based on 
the [Om]/H/3 and [Oll]/[Olll] ratios more uncertain, espe- 
cially for low O/H (see Fig. 10; we do not discuss here the 
[O ni]A4363 line, which is a diagnostic of low oxygen abun- 
dance but is very weak in most galaxies; see Kobulnicky et 
al. 1999 and Figs. 5-7 above). In case C, the unavailability 
of H/3 (and hence Ha/H/3) makes the uncertainties in e~ TV 
larger, but rjua and O/H can still be determined to within 
a factor of 2-3 using the [On], [Om], Ha, and [Nn] lines. 
As above, we find that the blending of Ha with the adjacent 
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Figure 10. Range in the logarithmic error A m (as defined by 
eq. fl9| ) obtained when applying the formulae in Table ^| to the 
models in Table |l| plotted against estimated line efficiency fac- 
tor, gas-phase oxygen abundance, and effective absorption optical 
depth of the dust. Each panel row corresponds to a different case 
in Table ^, as indicated. 



a galaxy are not expected to improve significantly SFR es- 
timates based on the [O n] luminosity. In nearby starburst 
galaxies, the correlation between ratio of far-infrared to ul- 
traviolet luminosities and ultraviolet spectral slope sets tight 
constraints on the effective absorption optical depth of the 
dust at ultraviolet wavelengths (e.g., Meurer, Heckman & 
Calzetti 1999; Chariot & Fall 2000). However, the different 
attenuation of line and continuum photons in these galax- 
ies (accounted for by values of greater than unity in 
Table |l]; see §2.4) make corrections of the [O il] luminos- 
ity uncertain. This is all the more critical in that the [O n] 
equivalent width is needed to estimate r]? n] (Table |^) . We 
therefore conclude that, without the Ha line, the SFR is 
difficult to estimate from the optical emission of a galaxy. 

It is worth recalling that the estimators in Table ^| ne- 
glect the potential contribution by active galactic nuclei to 
the ionizing radiation in galaxies. The presence of an AGN 
in a galaxy is usually revealed by readily identifiable sig- 
natures. These include larger [Oin]/H/3 and [Sn]/Ha ra- 
tios than in the galaxies in Fig. 2 and strong emission 
lines of highly-ionized species, such as [Nev]A3426 and 
[Fe vn]A6087, which are generally not produced by the softer 
spectra of stellar populations (e.g., Osterbrock 1989). Fur- 
thermore, the broad emission lines of most AGNs have ve- 
locity widths (FWHMs of several thousand kms~ ) that are 
several times larger than those expected from the virial mo- 
tions within galaxies (FWHMs less than about 1000 kms -1 , 
corresponding to line-of-sight velocity dispersions less than 
about 400 kms -1 ). Hence, the contamination of the nebular 
emission from a galaxy by an AGN can usually be identified. 



[N n] lines affects only weakly the accuracy in rjHa (case D). 
In the case where the H/3 line is observed but the [O n] line 
is not (case E), tjH„ and O/H can again be estimated to 
within a factor of 2-3 using the H/3, [Oin], Ha, and [Nil] 
lines, while the Ha/H/3 ratio helps constrain e~ rv . 

A critical situation arises when the Ha line is not avail- 
able (cases F and G in Table |^). We find that, in this case, 
the star formation rate cannot be determined reliably from 
the [On], H/3, and [Om] lines and the 4000 A discontinuity. 
Using all these features, we did not succeed in estimating 
??[On] ( or Vnp or J?[om]) to within less than a factor of 5 
(case F). In the case where only [On] and D4000 are avail- 
able (case G), the uncertainty in non] reaches a factor of 
almost 10, and O/H is unconstrained (Fig. 10). The major 
source of uncertainties in r?ron] estimates is the lack of a 
strong dust discriminant in the wavelength interval between 
[On] and [Om]. The uncertainties in rj? 0n , in cases F and 
G demonstrate that, without dust, the [On] efficiency fac- 
tor can be constrained to within a factor of only 2-3 based 
simply on the [O n] equivalent width and the 4000 A discon- 
tinuity. As a further check, we reestimated rjton] in case F 
after reducing the ranges in zero-age effective ionization pa- 
rameter and effective metallicity of the models in Table ^ 
to 0.5 Z e < Z < Z e and -3.0 < log((!7o)) < -2.0. The 
resulting uncertainty was 0.65, comparable to the value of 
0.70 in Table |^. We note in Fig. 10 that 7j[on] estimates are 
slightly more accurate for non] ^ 10 14 erg g _1 , correspond- 
ing to galaxies with high effective metallicity and high dust 
content, and ?7[oii] 10 15 erg g _1 , corresponding to galaxies 
with low effective metallicity and low dust content. 

Independent constraints on the absorption by dust in 



5 CONCLUSIONS 

We have developed a model for computing consistently the 
line and continuum emission from galaxies, based on an ide- 
alized description of the Hn regions and the diffuse gas 
ionized by single stellar generations in a galaxy. We have 
calibrated the nebular properties of this model using the 
observed [Om]/H/3, [O n]/[0 in], [Sn]/Ha, and [Nn]/[Sn] 
ratios of a representative sample of nearby spiral and irregu- 
lar, starburst, and Hn galaxies. To compute whole (line plus 
continuum) spectral energy distributions, we have included 
the absorption by dust in the neutral ISM using a recent 
simple prescription, which is consistent with observations of 
nearby starburst galaxies. The model succeeds in reproduc- 
ing quantitatively the optical spectra of nearby galaxies of 
various types. The prime parameters in our model, including 
the star formation history, the zero-age effective ionization 
parameter, the effective gas metallicity, the effective dust- 
to-heavy element ratio in the ionized gas, and the effective 
optical depth of the dust in the neutral ISM (and the fraction 
of this contributed by the ambient ISM) each have a specific 
influence on the integrated line and continuum properties of 
galaxies. Optical spectral fits enable us, in turn, to constrain 
simultaneously the star formation history, the metallicity, 
and the absorption by dust in galaxies. 

We have used this model to derive new estimators of 
the Ha and [On] efficiency factors (ijh and ??[oii]), the gas- 
phase oxygen abundance, and the effective absorption opti- 
cal depth of the dust in star-forming galaxies. The observed 
integrated spectral properties of nearby galaxies are com- 
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Table 2. Estimators of the Ha and [On] efficiency factors, the gas-phase oxygen abundance, and the effective absorption optical depth 
of the dust in galaxies for different assumptions about the available spectral information (cases A— G). 

Case Available spectral information^ Diagnostics^ and estimators^] 

[On] D 40 oo H/3 [Om] Ha [Nn] [Sn] 

A yes yes yes yes yes yes yes rjn a = 1.85 X f 15 x~ 2 ' 25 x^ 46 x% 90 X7 1 ' 14 X.T - 15 

O/H = 5.09 X 10- 4 x% 17 xl- 17 
e -rv = 1. 2n.Tr 3 - 32 .ri-06. r 0.62_ r -0.85_ r -0.49 



1.20 Xg Xg 2 X4 Xg 

^0 ^ = L44 x 10 15 3.0.75 x -0.17 x o.l5 x -0.58 
yes yes yes yes blended^] yes ma+[Nn] = 2.29 X 10 15 x\ s ° x\ 0A x~ 1A0 x 



O/H = 4.21 x 10- 4 x° 2 51 x 3 -°- 24 x 8 -°- 42 (L [0n] /L [0m] < 0.8) 
O/H = 4.58 x 10- 4 x-°' 37 x- a54 (L [0ll] /L [0m] > 0.8) 



2.01 „1.21 _-1.41 -1.70 
2 x 3 x 7 X S 



e~ T v = 0.90 x 

„0 _ n o, 1(1 15 „0.28 ..0.12 „-0.36 -0.06 

'Ha + [JVn] ~~ z - 01 A ±u J; 2 ^3 x 7 x 8 

yes yes no yes yes yes ... rin a = 1.36 X 10 15 xP/ 38 x ,- 80 x]~ ' 40 x ^ 33 

O/H = 3.98 X 10- 4 x" ' 21 x? 5 39 
e-*v =0.64xJ°' 12 4 62 x» 1 20 
r?« a = 1.69 x 10 15 x«- 64 x" 010 X" - 45 xO ™ 

yes yes no yes blended 6 ... r? Ha +[iV n] = 2.99 X 10 15 x~ 12 xj; 86 x° 3 68 

O/H = 4.32 Xl0- 4 x0- 13 x- - 43 (i[Oii]A[o m ] <0.8) 

O/H = 3.63 x 10- 4 x-°- 57 x-°- 72 (L [0ll] /L [0m] > 0.8) 

e-*v = 0.82x-°- 17 x° 12 61 x°- 3 35 

n° - 2 09 X TO 15 r° 14 r~° 13 

"Ha+[iVn] z - uy x iu x 13 x 14 

yes yes yes yes ... rm a = 1.71 X 10 15 x" 2 - 96 x~ 12 x~ 0A7 x" 18 X" - 38 

O/H = 4.15 x 10" 4 x" ' 29 x° 40 (i[Om]/^H3 < 2.5) 

O/H = 4.15 x 10- 4 X" °' 29 x% 40 (L[om]/Lng > 2.5) 

e -f v =1.15 x~ 3 ' 97 x~ 001 x" 1 ' 21 x ^ 01 x^f 2 
rj^ = 1.34 X 10 15 x 2 - 18 x" 012 x 2 - 44 x" 017 x" 2 ' 13 

yes yes yes yes r? [0 n] = 6.92 X 10 14 x~ °' 02 xg' 41 x^ 52 

O/H = 3.78 x 10 -4 Xj' 17 Xg -0,44 (i[On]A[Oni] <0.8) 

O/H = 3.96 xl0- 4 x-°- 46 (I[On]/i[Om] > 0.8) 

e-fv = .47x0- 05 x0- 01 

,» On] =1.43xl0^ 2 - 1M x-^x- ia x« 



G yes yes r? [0 n ] = 4.44 X 10 14 xfj 3 x\f 7 

O/H : 



e~ T v = 0.47 x^g 



-0.23 
17 
0.86 



„0 _i 91 „ ln 15 r 0.89 „3.78 
'[On] ~~ U 16 X 17 



Uncertaintyrl 


(logarithmic) 


0.26 


(0.39) 


0.24 


(0.26) 


0.27 


(0.46) 


0.14 


(0.20) 


0.33 


(0.39) 


0.52 


(0.58) 


0.21 


(0.26) 


0.32 


(0.38) 


0.25 


(0.26) 


0.41 


(0.52) 


0.31 


(0.40) 


0.49 


(0.54) 


0.12 


(0.15) 


u. oz 


l^U.OO ) 


0.52 


(0.56) 


0.27 


(0.41) 


0.50 


(0.59) 


0.25 


(0.30) 


U.40 


\\).0 ( ) 


0.32 


(0.33) 


0.40 


(0.40) 


0.35 


(0.49) 


0.21 


(0.35) 


0.70 


(0.99) 


0.51 


(0.57) 


0.31 


(0.31) 


0.55 


(0.56) 


0.20 


(0.40) 


0.94 


(1.47) 


0.63 


(0.63) 


0.56 


(0.61) 


0.39 


(0.96) 



a Observed line luminosities and D4000 (as defined by Bruzual 1983) not corrected for dust and stellar (H-Balmer) absorption. A 'yes' 
means that the observed rest-frame equivalent width of the line exceeds 5 A, a 'no' that the equivalent width is less than this limit, and 
an ellipsis that the line is not observed. The lines are [0 11] A3727, [O in] A5007, [N n] A6583 (see however footnote e), and [S 11] AA6717, 6731. 
6 The basis functions x^ are spectral features normalized to the same standard model as in cq. |l7) but including the effects of stellar 
Ha and H/3 absorption (for ry = 0): xi = (L Ha /iH/3)/3.4; x 2 = (L[o n]/ L [0 in] )/l-S; X3 = (L[ Om ]/L- a/3 )/2.0; x 4 = (L[ Sii ]/L Hq )/0.26; 
Xb = W Hj 9/(25A); x 6 = (L [Nii] /L [Sii] )/0.85; x 7 = (£H<*+[Nit]/-kH/3)/4.5; x g = (L [Sn] /L Ha+[Nn] )/0.20; x 9 = (L [Oii ]/L Hq )/0.90; x 10 = 
(i[Nn]/iHa)/0.22; xn = WW(170A); x 12 = (i[ O n]/^H«+[Nii])/0.r0; X13 = W Ha+[Nu] /(220 A); X14 = (L [Q In ] /L Ha +[N 11] )/0.45; 
xis = (^[onj/^H^J/S-O; xi6 = VK[ O ,i]/(60 A); xn = D 4 ooo/l-l- 

c Line efficiency factors jjhcu ^[Oii]! anc l *?Ha+[JVn] m units of erg g — 1 (eq. |16|); oxygen abundance by number relative to hydrogen in 
the ionized gas after accounting for the depletion onto dust grains [(0/H) cosm i c = 6.6 X 10 -4 ]; and effective dust-absorption factor e~ Tv 
at 5500 A (eqs. and 0). The line efficiency factors ify a , V?Q n y an d ^Ha+liVii] corres P on d to the case fy = 0. 

d 99th percentile range in |A m J (eq. ]^g|) obtained when applying the formula to all the models m in Table [l] (see text and Fig. 10). For 
reference, the full (100th percentile) range in |A m | is indicated in parentheses. 
e Ha+[Nn]AA6548, 6583 blend only (W Hq+ [jvii] > 5 A). 
© 1994 RAS, MNRAS 000, fj] [r] 
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patible with ranges of at least a decade in r/ua and Won]- 
The allowed combinations of ionized-gas parameters alone 
imply factors of 3.5 and 14 spreads in r/Ha and Won]) respec- 
tively. We find that, by exploiting the correlations between 
rjua and 7?[on] and various spectral features, we can reduce 
considerably the uncertainties affecting SFR estimates. For 
example, with the help of other lines such as [O n], H/3, [O ill] 
[Nil], or [Sn], the uncertainties in SFR estimates based on 
Ha can be reduced to a factor of only 2-3. A similar ac- 
curacy can be achieved even if the Ha line is blended with 
the adjacent [Nil] lines. Without Ha, however, the SFR, 
the gas-phase oxygen abundance, and the effective absorp- 
tion optical depth of the dust are difficult to estimate from 
the [O n] , H/3, and [O ill] lines. The reason for this is that 
the absorption by dust in the neutral ISM and the ionized- 
gas parameters are then difficult to constrain independently. 
This suggests that, while insufficient by itself, the Ha line 
is essential for estimating the star formation rate from the 
optical emission of a galaxy. 

The estimators listed in Table ^ allow one to derive 
more refined constraints on the star formation rate, the gas- 
phase oxygen abundance, and the effective absorption opti- 
cal depth of the dust in galaxies than was previously possi- 
ble from optical spectra alone. By design, these estimators 
can be applied directly to observations uncorrected for dust 
and stellar (H-Balmer) absorption. We can also derive, on 
demand, the estimators best suited to observational selec- 
tion criteria other than those exemplified in Table ^. Finally, 
as shown in §3.2, our full model for the line and continuum 
emission allows detailed interpretations of whole spectral en- 
ergy distributions of galaxies in terms of stars, gas, and dust 
parameters. 



ACKNOWLEDGMENTS 

We thank M. Fall, G. Ferland, R. Kennicutt, and L. Tresse 
for valuable discussions. S.C. appreciates the generous hos- 
pitality of the MPIA-Garching, and M.L. acknowledges 
support by the European Commission under TMR grant 
no. ERBFMBI-CT97-2804. This research was supported in 
part by the National Science Foundation through grant 
no. PHY94-07194 to the Institute for Theoretical Physics. 

REFERENCES 

Allard F., Hauschildt P. H., 1995, ApJ, 445, 433 

Alongi M., Bertelli G., Bressan A., Chiosi C, Fagotto F., Greggio 

L., Nasi E., 1993, A&AS, 97, 851 
Baldwin J. A., Phillips M. M., Terlevich R., 1981, PASP, 93, 5 
Barbaro G., Poggianti B. M., 1997, A&A, 324, 490 
Bessell M. S., Brett J. M., Scholtz M., Wood P. R., 1989, A&AS, 

77, 1 

Bessell M. S., Brett J. M., Scholtz M., Wood P. R., 1991, A&AS, 
89, 335 

Binette L., Magris C. G., Stasinska G., Bruzual A. G., 1994, A&A, 
292, 13 

Bottorff M., LaMothe J., Momjian E., Verner E., Vinkovic D., 

Ferland G., 1998, PASP, 110, 1040 
Bresolin F., Kennicutt R. C, Garnett D. R., 1999, ApJ, 510, 104 
Bressan A., Fagotto F., Bertelli G., Chiosi C, 1993, A&AS, 100, 

647 

Bressan A., Chiosi C, Tantalo R., 1996, A&A, 311, 425 



Bruzual A. G., 1983, ApJ, 273, 105 

Bruzual A. G., Chariot S., 1993, ApJ, 405, 538 

Calzetti D., 1997, AJ, 113, 162 

Calzetti D., Kinney A. L., Storchi-Bergmann T., 1994, ApJ, 429, 
582 

Campbell A., 1988, ApJ, 335, 644 

Cananzi K., Augarde R., Lequeux J., 1993, A&AS, 101, 599 

Cassinelli J. P., et al., 1995, ApJ, 438, 932 

Chariot S., Fall S. M., 1993, ApJ, 415, 580 

Chariot S., Fall S. M., 2000, ApJ, 539, 718 

Cowie L. L., Songaila A., 1986, ARAA, 24, 499 

Diaz, A. I., Terlevich, E., Vflchez, J. M., Pagel, B. E. J., Edmunds, 

M. G., 1991, MNRAS, 253, 245 
Draine B. T., Lee H. M., 1984, ApJ, 285, 89 
Evans I. N., Dopita M. A., 1985, ApJS, 58, 125 
Fagotto F., Bressan A., Bertelli G., Chiosi C, 1994a, A&AS, 100, 

647 

Fagotto F., Bressan A., Bertelli G., Chiosi C, 1994b, A&AS, 104, 
365 

Fagotto F., Bressan A., Bertelli G., Chiosi C, 1994c, A&AS, 105, 
29 

Fanelli M. N., O'Connell R. W., Thuan T. X., 1988, ApJ, 334, 
665 

Ferguson A. M. N., Wyse R. F. G, Gallagher J. S., Hunter D. A., 

1996, AJ 111, 2265 
Ferland G., 1996, Hazy, A Brief Introduction to Cloudy. Internal 

Report, Univ. Kentucky, USA 
Fluks M. A., Plez B., The P. S., de Winter D., Westerlund B. E., 

Steenman H.C., 1994, A&AS, 105, 311 
Garnett D. R., Dufour R. J., Peimbert M., Torres-Peimbert S., 

Shields G. A., Skillman E. D., Terlevich E., Terlevich, R. J., 

1995, ApJ, 449L, 77 
Girardi L., Bressan A., Chiosi C, Bertelli G., Nasi E., 1996, 

A&AS, 117, 113 
Gallagher J. S., Hunter D. A., Bushouse H, 1989, AJ, 97, 700 
Garcia- Vargas M. L., Bressan A., Diaz A. I., 1995, A&AS, 112, 

35 

Giallongo E., Fontana A., Madau P., 1997, MNRAS, 289, 629 
Gonzalez-Delgado, R.M., Leitherer, C, 1999, ApJS, 125, 479 
Grevesse N., Anders E., 1989, in Waddington C.J., ed., Proc AIP 

Conference 183, Cosmic abundances of matter. New York, p. 

1 

Grevesse N., Noel A., 1993, in Prantzos N., Vangioni-flam E., 
Casse M., eds, Origin & Evolution of the Elements. Cam- 
bridge Univ. Press, p. 15 
Henry R. B. C, Worthey G., 1999, PASP, 111, 919 
Hummer D. G., Mihalas D., 1970, MNRAS, 147, 339 
Hummer D. G., Storey P., 1992, MNRAS, 254, 277 
Hoopes C. G., Walterbos R. A. M., Greenwalt B. E., 1996, AJ, 
112, 1429 

Hunter D. A., Gallagher J. S, 1990, ApJ, 362, 480 
Izotov Y. I., Thuan T. X., Lipovetzky V. A., 1994, ApJ, 435, 647 
Jenkins E. B., 1987, in Hollcnbach D., Thronson H, eds, Inter- 
stellar Processes. Reidel, Dordrecht, p. 533 
Kennicutt R. C, 1983, ApJ, 272, 54 
Kennicutt R. C, 1984, ApJ, 287, 116 
Kennicutt R. C, 1992a, ApJS, 79, 255 
Kennicutt R. C, 1992b, ApJ, 388, 310 
Kennicutt R. C, 1998, ARA&A, 36, 189 

Kennicutt R. C, Tamblyn P., Congdon C. E., 1994, ApJ, 435, 22 
Kinney A. L., Bohlin R. C, Calzetti D., Panagia N., Wyse R. F. 

G., 1993, ApJS, 86, 5 
Kobulnicky H. A., Kennicutt R.C., Pizagno J.L. 1999, ApJ, 514, 

544 

Kurucz R. L., 1992, in Barbuy B., Renzini A., eds, Proc. IAU 
Symp. 149, The Stellar Populations of Galaxies. Kluwer, Dor- 
drecht, p. 225 

Kurucz R. L., 1995, private communication to T. Lejeune 



© 1994 RAS, MNRAS 000, 



Nebular Emission from Star-Forming Galaxies 17 



Leitherer C, Heckman T. M., 1995, ApJS, 96, 9 

Leitherer C, Ferguson H. C, Heckman T. M., Lowenthal J. D. 

1995, ApJ, 454L, 19 
Lehnert M. D., Heckman T. M., 1994, ApJ, 426, 27 
Lejeune T., Cuisinier F., Buser R., 1997, A&AS, 125, 229 
Lejeune T., Cuisinier F., Buser R., 1998, A&AS, 130, 65 
Martin C. L., 1997, ApJ, 491, 561 

Martin P. G., Rouleau F., 1989, in Malina R. F., Bowyer S., eds, 
Extreme Ultraviolet Astronomy. Pergamon, Elmsford, NY, p. 
34 

McCall M. L., Rybski P. M., Shields G. A., 1985, ApJS, 57, 1 
McGaugh S. S., 1991, ApJ, 380, 140 

McQuade K., Calzetti D., Kinney A. L., 1995, ApJS, 97, 331 
Mathis S. J., 1986a, PASP, 98, 995 
Mathis S. J., 1986b, ApJ, 301, 423 

Meurer G. R., Heckman T. M., Calzetti D., 1999, ApJ, 521, 64 
Oey M. S., Kennicutt R. C, 1993, ApJ, 411, 137 
Oey M. S., Kennicutt R. C, 1997, MNRAS, 291, 827 
Osterbrock D. E., 1989, Astrophysics of Gaseous Nebulae and 
Active Galactic Nuclei. Mill Valley: University Science Books 
Scalo, J.M., 1986, FCPh, 11, 1 
Schaerer D., de Koter A., 1997, A&A, 322, 598 
Shields J. C, Kennicutt R. C, 1995, ApJ, 454, 807 
Shields G. A., Searle L., 1978, ApJ, 222, 821 
Snow T. P., Witt A. N., 1996, ApJ, 468L, 65 

Spitzer L., 1978, Physical Processes in the Interstellar Medium. 

Wiley, New York, p. 333 
Stasirlska G, 1980, A&A, 85, 359 
Stasinska G, 1990, A&AS, 83, 501 
Stasinska G, Leitherer C, 1996, ApJS, 107, 661 
Storchi-Bergmann T., Kinney A.L., Challis P., 1995, ApJS, 98, 

103 

Terlevich R., Melnick J., Masegosa M., Moles M., Copetti V.F., 

1991, A&AS, 91, 285 
Tresse L., Rola C, Hammer F., Stasinska G., Le Fevre O., Lilly 

S. J., Crampton D., 1996, MNRAS, 281,847 
Vacca W. D., 1994, ApJ, 421, 140 

Vacca W. D., Garmany C. D., Shull J. M., 1996, ApJ, 460, 914 
Vila-Costas, M. B., Edmunds M. G, 1993, MNRAS, 265, 199 
Wang J., Heckman T. M., Lehnert M. D., 1997, ApJ, 491, 114 
Worthey G., Faber S. M., Gonzalez J. J., BursteinD., 1994, ApJS, 
94, 687 

Zaritsky D., Kennicutt R. C, Huchra J. P., 1994, ApJ, 420, 87 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science I^TjtX style file. 



© 1994 RAS, MNRAS 000, §-0 



